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Abstract 

The possibility of measuring neutral-current coherent elastic neutrino-nucleus 
scattering (CENNS) at the TEXONO experiment has opened high expectations to¬ 
wards probing exotic neutrino properties. Eocusing on low threshold Germanium- 
based targets with kg-scale mass, we find a remarkable efficiency not only for de¬ 
tecting CENNS events due to the weak interaction, but also for probing novel 
electromagnetic neutrino interactions. Specifically, we demonstrate that such ex¬ 
periments are complementary in performing precision Standard Model tests as well 
as in shedding light on sub-leading effects due to neutrino magnetic moment and 
neutrino charge radius. This work employs realistic nuclear structure calculations 
based on the quasi-particle random phase approximation (QRPA) and takes into 
consideration the crucial quenching effect corrections. Such a treatment, in con¬ 
junction with a simple statistical analysis, shows that the attainable sensitivities 
are improved by one order of magnitude as compared to previous studies. 
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1 Introduction 


The robust discovery of neutrino oscillations in the propagation of solar and atmospheric 
neutrinos BE], conhrmed at accelerator and reactor neutrino sources PS] has provided 
us with a rather solid proof for the existence of neutrino masses and mixing EE] and 
hence the clearest evidence for the need of physics beyond the Standard Model (SM) [7]. 
These results have prompted a great rush to produce adequate SM extensions with small 
neutrino masses |8]. 

Underpinning the ultimate origin of neutrino mass stands out as one of the biggest 
challenges in particle physics [H]. A generic feature of many such schemes is the presence 
of non-vanishing neutrino electromagnetic (EM) properties [IHl HH [I2l [ISl UTl US] . While 
the neutrino masses indicated by oscillation data are perhaps too small to induce sizeable 
magnetic moments, this issue is rather model dependent, and one cannot exclude this 
possibility on general grounds [16]. If large enough, these may still play an important sub¬ 
leading role in precision neutrino studies nil, despite the good agreement found within 
the three-neutrino oscillation picture. Non-zero diagonal magnetic moments exist for 
massive Dirac neutrinos. In contrast, in the general Majorana neutrino case all magnetic 
moments are transition-type. Therefore, the study of neutrino magnetic moments would 
be a powerful tool towards distinguishing their Dirac or Majorana character [T8l IT9] . 

The detection of neutral-current (NC) coherent elastic neutrino-nucleus scattering 
(CENNS) processes by measuring the nuclear recoil spectrum of the scattered nucleus 
has by now become feasible IZD]. As a concrete example, the newly formed COHERENT 
Collaboration at the Spallation Neutron Source (SNS) has excellent prospects [211 [22], 
motivating also theoretical effort [23l [23]. In this work, we consider the possibility of 
revealing signs of new physics through a detailed study of CENNS at the TEXONO 
experiment [231 ESI EZ]. We demonstrate that the use of sub-keV Germanium-based 
kg-scale detectors [2EI EH] ED], provides a favourable experimental set up with good 
prospects for performing precision SM tests, as well as probing EM neutrino proper¬ 
ties [31] , such as the neutrino magnetic moment [TH [191 EH [33] and the neutrino charge 
radius [331 ESI ESI EH EH] . The total number of events expected in an experiment search¬ 
ing for CENNS depends strongly on the energy threshold Tth^es, as well as the total mass 
of the detector [231123]. For low energy thresholds and more massive detectors, the total 
number of events expected is significantly larger and, therefore, the attainable sensitiv¬ 
ities are higher. In addition, this work highlights that the present calculations become 
more realistic by considering quenching effect corrections [39l [301 [31] . The sensitivity is 
evaluated by assuming that a given experiment searching for CENNS events will measure 
exactly the SM expectation. Thus, any deviation [3H SEl SI] is understood as a signature 
of new physics [3S1 SS] SH SSI SH] . 
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Apart from the possibility of the hrst ever detection of CENNS events, our present 
results emphasise the potentiality of discovering neutrino interactions beyond the SM 
expectations [50]. In our estimates we perform nuclear structure calculations within the 
context of the quasi-particle random phase approximation (QRPA) that uses realistic nu¬ 
clear forces and employ a x^-type statistical analysis. We hnd that the 

prospects for improving current bounds on are rather promising and complementary 
to future sensitivities on the muon neutrino magnetic moment, m 

2 Coherent elastic neutrino-nucleus scattering 

The coherent elastic scattering of neutrinos upon a nucleus is described within the SM 
starting from the neutrino-quark NC interaction Lagrangian. However, as mentioned 
above, one has good reasons to expect corrections coming from new physics [I], such as 
non-standard interactions (NSI) |l2l HS] HU US] HHl HU HQ] or non-trivial neutrino EM 

properties [ini[ni[I21[I31[Il[il[ISlE2lE31ES|. 

2.1 Standard model prediction 

Within the context of SM, for low energies {Ey -C My/) accessible to neutrino exper¬ 
iments, the weak neutral-current CENNS can be naturally studied by considering the 
V E A interaction of four-fermion vvff type operators entering the effective Lagrangian, 
written as [13| 


/:SM = -2V2Gp Y. si’a [i^alpLUa] [f^Pf] • (1) 

P=L,R 
f= u^d 
a.=e,fi^T 

In this Lagrangian, the chiral structure of the SM weak interaction is expressed by the 
projectors P = {L,R}, while a denotes the neutrino flavour and / refers to the hrst- 
generation quarks. The relative strength of the left- and right-handed couplings for u- 
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and d-quark to the Z-boson with respect to the Fermi constant Gp are given as 


_„NC / 1 _ 2 - a2 1 _L \'^,L 
9 act —PuN I 2 ^ ’ 

_ NC ( ^ I , s2 \ , td,L 

9aa —PuN I “2 ’ 

^u,R _NC f 2 -2 \ I \u,R 

9aa —PuN \~^I^vNSz ) + ’ 

^d,R_Nc(^- -2\ I \d,R 

9aa —PvN 1 -^'^i'NSz ) + 


(2) 


In the latter expressions, after inclnding the relevant radiative corrections, we have s| = 
sin^^vK = 0.23120, = 1.0086, k^N = 0.9978, = -0.0031, A'^’^ = -0.0025 and 

^d,R _ 2X^’R = 7.5 X 10“^ [IT]. From the effective Lagrangian in Eq. Q, the differential 
cross section with respect to the nnclear recoil-energy, T, for the case of a CENNS off a 
spherical spin-zero nnclens of mass M, reads 


da A 

'^^7 SM 


G%M 

271 


MT 


. T 


(^s||7Wo(g)||^s) 


(3) 


For gs —)■ gs transitions, the corresponding coherent nnclear matrix element takes the 
form jlHl [23| 


{gs\\Mo{q)\\gs) = [2{g^’^ -f g^’^) + {g'^^ + g'^^)] ZFz{q^) 

+ [{Pat + 9aa) + + S'aif)] NF^i^q^) ■ 

Note that, dne to the smallness of the conpling of protons with the Z-boson, the main 
contribntion to the CENNS cross section essentially scales with the sqnare of the nentron 
nnmber N of the target nnclens (see e.g. [M]). We stress that the differential cross 
section is evalnated with high signihcance by weighting the nnclear matrix element with 
corrections provided by the proton (nentron) nnclear form factors Fz(N){.q^)- This way the 
hnite nnclear size is taken into acconnt with respect to the typical momentnm transfer, 
q ~ a/2MT. Furthermore, the enhancement of the CENNS cross section makes the 
relevant experiments favourable facilities to probe the neutron form factor of the target 
nucleus at low energies [20l|2ll[22|. 

From a nuclear theory point of view, the reliability of the present CENNS cross 
sections calculations is maximised in terms of accuracy by performing nuclear structure 
calculations in the context of QRPA [52[ |53]. Motivated by its successful application 
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on similar calculations for various semi-leptonic nuclear processes (SH ES], in this work 
we construct explicitly the nuclear ground state, \gs) = |0'''), of the relevant even-even 
isotope through the solution of the BCS equations (for a detailed description see Ref. [23]). 


2.2 Electromagnetic neutrino-nucleus cross sections 

The existence of neutrino masses is well-established thanks to the current neutrino oscil¬ 
lation data, implying that they could have exotic properties, such as non-zero neutrino 
magnetic moments. In this framework, potential neutrino-nucleus interactions of EM 
nature have been considered resulting in corrections to the weak 

CENNS cross section of the form [3l] 

(^) -(%) +(%) . (5) 

Here, the helicity-violating EM contribution to the neutrino-nucleus cross section can be 
parametrised in terms of the proton nuclear form factor, the hne structure constant Oem 
and the electron mass rUg as [I6] 

Note, that additional corrections are incorporated in Eqs. (§, (§, compared to Ref. [21] . 

If the neutrino is of Dirac-type as in the SM, then the magnetic moment is un- 
detectably small due to its proportionality to the neutrino mass. Even though from 
oscillation data the latter is well-known to be small, one cannot rule out the possibility 
of sizeable neutrino magnetic moments. Indeed, in general scenarios where neutrinos are 
Majorana fermions, as expected on general grounds, larger transition magnetic moments 
are possible. For example, relatively sizeable contributions may be predicted in models 
involving NSI [16]. From experimental perspectives, a potential signal will be detected as 
a distortion of the nuclear recoil spectrum at very low energies where the EM cross section 
dominates due to its ~ 1/T dependence. For this reason, such challenging technological 
constraints require innovative experimental advances towards reducing the threshold to 
the sub-keV region. 

Despite having vanishing electric charge, the hrst derivative in the expansion of the 
neutrino electric form factor entering the decomposition of the leptonic matrix element 
may provide non-trivial information concerning other neutrino electric properties [15] . 
More specihcally, the gauge-invariant dehnition of the neutrino root mean square charge- 
radius a = {e, p, r} [3ll|35] leads to corrections to the weak mixing angle [36llT7] l38] 
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3 The TEXONO experiment 


In the present work we explore how well one can probe neutrino EM phenomena with the 
TEXONO experiment |25l |26l [27] through low-energy CENNS measurements near the 
Kuo-Sheng Nuclear Power Station. Towards this purpose, the TEXONO Collaboration 
has pursued a research program aiming at detecting neutrino-nucleus events by using high 
purity Germanium-based detectors HPGe with sub-keV threshold [2Sl[2Hll3n]- According 
to the proposal, we consider a 1 kg ^®Ge-detector operating with a threshold as low 
as Tthres = 100 eVee- Due to the absence of precise information regarding the fuel 
composition of the reactor core, we only include the dominant component of the anti¬ 
neutrino spectrum. In this respect, for the present study we assume a typical neutrino 
flux of = 10^^ for a detector location at 28 m from the reactor core. In 

order to estimate the reactor antineutrino energy-distribution rip^{E^) for energies above 
2 MeV, existing experimental data from Ref. [56] are employed. We stress that the 
main part of reactor anti-neutrinos is released with energies < 2 MeV, thus their 
contribution is crucial and must be taken into account. For their description we adopt 
the theoretical estimates given in Ref. |S7|. This will bring about improved sensitivities 
on the neutrino magnetic moment. 

4 Results and discussion 

4.1 Signal cross sections 

First we present and discuss the individual weak and electromagnetic differential and 
total GENNS cross sections (see Eq. weighted over experimental reactor anti-neutrino 
spectra [SniET]. These convoluted cross sections determine the neutrino signals expected 
to be recorded at a nuclear detector (e.g. the ^®Ge of the TEXONO experiment). For 
each interaction channel x, [x = SM, EM, tot] the energy-integrated differential cross 
section, {da/dT)^, is dehned as 
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Figure 1: Weak and electromagnetic differential {left panel) and total {right panel) cross 
sections convoluted with reactor z/g-spectra. 


where r]p^{Ei,) denotes the normalised neutrino energy-distribution. The corresponding 
signal cross section reads 

<»■>(£„) = jdT{^(E,,T)) (9) 

For the ^®Ge detector, assuming the experimental constraints placed recently by TEX- 
ONO (/ip^ = 7.4 X 10“^^/is [25] and = 6.6 x 10“^^cm^ |26|), the computed results 
are illustrated in Fig. One sees that in the case of the curve involving neutrino 
magnetic moment contribution exceeds that of the pure SM weak rate at low neutrino 
energies, Ey. The curve containing neutrino charge radius contributions through Eq. 
is showing a similar behaviour as the pure SM (the photon propagator cancellation leads 
to 4-fermion contact interaction [30l l50] ). 

4.2 Statistical analysis 

Our present analysis is strongly based on the estimation of the number of CENNS events. 
Therefore, we hrst provide a brief description of the conventions and approximations we 
use in our calculations. For each interaction channel x, the number of CENNS events 
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Figure 2: sensitivity profile as a function of the weak mixing angle 9\y at TEX- 

ONO. The results are presented for two detector efficiencies 100% (50%) shown with red 
(black) colour and three values of the quenching factors Qj = (1,0.20,0.25) indicated 
with (dashed, solid, dashdotted) lines. (For interpretation of the references to colour in 
this hgure legend, the reader is referred to the web version of this article.) 


above a minimum nuclear recoil-energy, Tmin, reads 


N, = K 


r]pXEu)dE^ 


%(K,T)] iT. 


( 10 ) 


In the above expression, K = Atargttot'hp^, with Atarg being the total number of atoms 
in the detector and ftot the relevant irradiation period. Note that potential effects due 
to neutrino oscillation in propagation are neglected, since this is well satished for the 
short-baselines considered here. The numerical results throughout this work refer to a 
1 kg ^®Ge-detector, one year of data taking and a detector threshold of 100 eVge. In 
addition, we consider two different detector efficiencies including an optimistic approach 
of a perfectly efficient detection capability and the more realistic scenario assuming a 
recoil acceptance of 50%. 

The present calculations take into consideration the fact that the nuclear recoil events 
are quenched [39] (in Ref. m, where for each target the calculation is referred to the 
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Table 1: Summary of the sensitivities obtained for sin^ 6w (la) and for the EM neutrino 
parameters (90% C.L.) at the TEXONO experiment. The results refer to various sen¬ 
sitivities and quenching factors. Comparing with Ref. one sees that a substantial 
improvement in the sensitivity for the weak mixing angle the magnetic moment 


parameter and the neutrino charge-radius (r?^) w.r.t. the COHERENT proposal. 


(Target, Threshold) 

EfBciency 

Quenching 

COHERENT’ 

(100 kg ’■"Go, 10 keVee 

) 


TEXONO (this work) 

(1 kg ^'‘Ge, 100 eVee) 



67% 

Qf = 1 


100% 



50% 


Q/ = 1 

Qf = 0.20 

Qf = 0.25 

Q/ = l 

Qf = 0.20 

Qf = 0.25 


0.0055 

0.0010 

0.0033 

0.0025 

0.0014 

0.0046 

0.0035 

Uncer. (100%) 

2.36 

0.43 

1.41 

1.08 

0.61 

1.97 

1.51 

X 10“^“ flB 

9.46 

0.40 

0.98 

0.83 

0.47 

1.17 

0.99 

(r?) X 10“®^ cm^ 

-0.38 - 0.37 

-0.07-0.07 

-0.22 - 0.22 

-0.17-0.17 

-0.10 - 0.10 

-0.32 - 0.31 

-0.24 - 0.24 


nuclear recoil-energy window [ 21 ], such a treatment is not necessary). Corrections of this 
type are crucial since for a given ionisation detector the observed energy (equivalent to an 
electron energy) is lower than the total nuclear recoil-energy, i.e. much energy is converted 
to heat (phonons) which is not measured, especially at low energies |10|. To convert from 
nuclear recoil-energy (eVnr) to electron equivalent energy (eVee), we multiply the energy 
scale by a quenching factor, Qf. In principle Qf varies with the nuclear-recoil energy |TT| 
and has to be determined experimentally, however for the sub-keV Germanium-based 
targets considered here it can be well-approximated as constant with typical values in 
the range 0.20-0.25 [28]. Thus, the TEXONO threshold Tthres = 100 eVee corresponds to 
nuclear recoil-energy T-^ain = 500 eVnr ioi Qf = 0.20 and Tmin = 400 eVnr for Q/ = 0.25, 
correspondingly the maximum nuclear recoil-energy T^ax = 1-81 keVnr, is restricted to a 
maximum observable energy of 362 and 452 eVee- 

The sensitivity of the TEXONO experiment to the SM weak mixing angle sin^ 9w is 
quantitatively determined on the basis of a y^-type analysis relying on statistical errors 
only [21] 


y 


2 


/ Nsm — -^SM(sin^ Oyi/) 
V SNsm 


2 


( 11 ) 


In our calculational procedure we have assumed that the TEXONO experiment will 
detect the precise number of SM events, Nsu, by hxing the electroweak mixing parameter 
to the PDG value i.e. s| = 0.23120 [I?]. We estimate the expected events Nsm = 
(27962, 2586, 4415) assuming Qf = (1, 0.20, 0.25) and a detection threshold 100 eVee, 
in good agreement with Ref. [29]. The y^ function is then minimised with respect to 
siri^ 9w, by varying this parameter around its central value, taken as the value reported 
by the PDG. 
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Figure 3: {Left panel) Ax^ profiles in terms of the neutrino magnetic moment parameter 
in units of at TEXONO. {Right panel) Sensitivity to /ip^ at 90% C.L. as a 

function of the detector mass. Same conventions as in Fig. [^are used. 


We have explicitly verihed that there are good prospects for making precision tests 
of the SM by using low-energy '^®Ge detectors. Our results for the TEXONO sensitivity 
to the weak mixing angle are presented in Fig. Furthermore, we have also evaluated 
the Icr error band on sin^ 6 ^ 1 ^ dehned as (5sin^6'vi/ = = (s^max — s‘^^i„)/2 as well as 

the corresponding uncertainty with s'^max being the respective upper 

(lower) bound. The resulting sensitivities are shown in Table Specihcally, neglecting 
the quenching corrections {Qf = 1), the improvement upon previous results [ 2 l] is up 
to 82% (74% when realistic efficiencies are taken into account). Furthermore, the effect 
ioT Qf = 0.20 (0.25) leads to reduction of 5sin^ 9iy sensitivity by a factor of 3.3 (2.5) for 
both detection efficiencies. 

Prompted by the upcoming generation of low-threshold nuclear detectors, we have 
made an effort to identify possible deviations from the SM neutrino-quark interaction 
cross section originated by non-standard neutrino EM properties. In particular, analysing 
their sensitivity to electromagnetic CENNS events we have found that important devia¬ 
tions may be induced by the presence of a non-zero transition neutrino magnetic moment 
fieff = Pi/e- In order to determine this sensitivity we use a function of the form [19] 

^ ( 12 ) 

V oNsm J 
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Figure 4: 90% C.L. allowed regions in the {siv? Ow-IJ^u^) plane {left panel) and the {{rlf}- 
fipj plane {right panel) from a two parameter combined analysis. See the text for more 
details. 


In Eq. (12), we substitute the SM cross section by the one given in Eq. (|^ in order to 
account for possible events, originating from the corrections associated to the non¬ 
trivial structure of the neutrino EM current, as discussed previously. The corresponding 
results obtained by varying the effective transition neutrino magnetic moment, /ip^, are 
presented in Fig. (left panel). The experimental TEXONO limit from he — e scattering, 
= 7.4 X [25], is also shown for comparison (the most stringent bound on 

the neutrino magnetic moment comes from the reactor experiment GEMMA as = 
2.9 X at 90% C.L. [HH])- One sees that the prospects are very promising. Indeed, 

from Table [T| we find that the attainable sensitivities are improved by about one order 
of magnitude compared to the corresponding expectations at a SNS facility, considered 
recently in [2l]. Note however, that experiments at the SNS are not optimised to measure 
electron-neutrino properties. 

Moreover, it is worth mentioning that for the case of a 100 eVee threshold and Qf = 
1, the resulting sensitivity is by 46% (36% for the case of realistic efficiency) better 
than the existing limits derived from he — e scattering TEXONO data 123. However, 
assuming Qj = 0.20 (0.25) the above sensitivity reduces by a factor of 2.5 (2.1) for both 
recoil acceptances. From our calculations we have also found that neglecting quenching 
corrections the sensitivity to /ip^ of a given detector with mass m is roughly equivalent 
to that of a detector with ten times bigger mass for the case of Qf = 0.25. The results 
concerning this point are shown in Fig. (right panel). 


10 




























In view of our previous discussion, the TEXONO sensitivity to {rp^)-related searches 
is estimated through the dehnition of the given in Eq. ( [TT| ) by replacing sin^ 9w with 
that of Eq. Q and hxing sin^ 9w to the PDG value. After the minimisation we 
hnd that the TEXONO experiment is expected to be very sensitive to EM contributions 
of this type. The estimated 90% C.L. sensitivities are presented in Table In the 
particular case of Qf = 0.25 we see that thanks to the observation of CENNS events, 
TEXONO can reach an improvement of the order of 35% or more, with respect to similar 
calculations m- Again, this sensitivity reduces by a factor of 3.2 (2.4) when quenching 
corrections Qf = 0.20 (0.25) are taken into account. Moreover, it is worth noting that 
the latter sensitivities are by one order of magnitude better than the current TEXONO 
constraint obtained from — e. scattering [22] . 

Finally, it is also interesting to show the combined sensitivities obtained by varying two 
of the above parameters (sin^^vi/, and (r?J) simultaneously. The 90% C.L. allowed 
regions in the (sin^ and {{rlQ-fipJ plane are shown in the left and right panel 

of Fig. 1^ respectively for different quenching factors. One notices that the resulting 
parameter space is substantially reduced with respect to the corresponding sensitivity 
regions for muon neutrinos at a SNS experiment (see e.g. Ref. 1231). The latter is 
a direct consequence of the low-threshold TEXONO detectors adopted in the present 
study. 

5 Summary and conclusions 

In this work we have explored the possibility of performing Standard Model precision 
studies and probing for new physics through low energy neutral-current coherent elastic 
neutrino-nucleus scattering (CENNS) at the TEXONO experiment. Moreover, we have 
presented a comprehensive analysis for the case of potential sub-leading neutrino EM 
interactions. The calculated convoluted cross sections, clearly indicate the need for novel 
detector technologies with sub-keV sensitivities. Furthermore, from a nuclear physics 
point of view, the reactor neutrino beam induces transitions in the bound nuclear spec¬ 
trum while the Spallation Neutron Source (SNS) beam may excite in addition much 
higher transitions of the nuclear detector. We conclude that, apart from providing the 
hrst ever detection of CENNS events, low threshold Cermanium-based kg-scale detec¬ 
tors, e.g. TEXONO, will bring substantial improvements on precision SM tests as well 
as sensitivities on neutrino EM properties, such as the neutrino magnetic moment and 
the neutrino charge radius. We show explicitly that the sensitivities improve by up to 
one order of magnitude with respect to previous estimates. In this paper we have used 
realistic nuclear structure calculations within the context of the quasi-particle random 
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phase approximation (QRPA) and a simple analysis taking into account quench¬ 

ing effects. We have also checked that our sensitivities are determined mainly by the 
number of events: a binned sample would result in differences less than percent. It is 
also worth mentioning that a global fit including an experiment at the SNS added to the 
TEXONO experiment will not essentially improve these results, since SNS provides the 
best sensitivities for the case of muonic neutrinos. Hence, the two experiments are clearly 
complementary. 
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